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Abstract 

We discuss the Skyrme model for strong interactions and the concept of a Skyrmion 
^SJ ' fluid. The corresponding compact objects, namely Skyrmion stars, constructed from 

^ ■ the equation of state describing such a fluid are presented and compared to models 

, of neutron stars based on modern equation of states. We suggest plausible Skyrmion 

I star candidates in the 4U 1636—53 and 4U 1820—30 low mass X-ray binary sys- 

^ ■ tems where the suggested masses of the accreting compact companion {^^ 2.0Mq) 

, remain a challenge for neutron star models. 

2 : 1 Introduction 

ci • We are closing in on neutron stars both observationally and theoretically. Observationally 
more and better masses and radii are determined by a number of different methods. The 
^ ! "large" masses implied in few cases (as in the Vela X-1 pulsar where a mass of ~ l.TSM© 

^ I has been suggested; poses a real challenge to models of neutron stars built using the 
so-called modern equation of state (EOS) where the uncertainties are reduced by improved 
two- and three-body forces, relativistic effects and many-body calculations |2j. Even the 
stiffest EOS so far developed seem to be facing difficulty in accounting for the extreme 
values (up to 2.4:Mq; Sect. 4). In [Hj, hereafter OB, as an alternative, we constructed 
an EOS of dense matter based on the Skyrme model for strong interactions which repre- 
sents baryons as solitons of classic pionic fields. The resulting compact objects we named 
Skyrmion stars (SSs) ^ are intrinsically heavier (due to the stiffness of the Skyrmion fluid; 
hereafter SF) than any other type of compact stars computed using modern EOS. SSs can 
be as massive as ~ 2.8Mq leading us to speculate (given the above mentioned observations) 
that these might exist in nature^. The paper is presented as follows: Sect. 2 is devoted 
to the discussion of the Skyrme model and its link to Quantum-ChromoDynamics (QCD). 

^Email: ouyed@phas.ucalgary.ca 

^Skyrmion stars, like neutron stars, are likely to be born with masses around 1.5Mq. We expect only 
older SSs that have accreted enough mass to reach these extreme masses. 
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The problem of the missing attraction in the Skyrme model is described. In Sect. 3, we 
discuss the role the dilaton (the glueball potential in QCD) could play in curing such a 
problem by binding Skyrmions together to form the SF. We end Sect. 3 by reminding the 
reader of the basic properties of the SF and the resulting EOS used to construct models 
of SSs. The astrophysical implications follow in Sect. 4 where we compare SSs to stars 
constructed using modern EOS of dense matter. We conclude in Sect. 5. 

2 The Skyrme model 

We first give a brief overview of the Skyrme model. The fundamental principles are dis- 
cussed at a basic level. The interested reader is referred to Ref. |2i for a thorough intro- 
duction to the topic. 

Skyrme constructed a model of pion interactions consisting of a conventional model 
of weak meson interactions plus an additional (higher-order) term thought to take into 
account indirect effects of heavier mesons like the p-meson. The now well-known Skyrme 
Lagrangian density is usually written as [U], 

-^Skyrme = L2 + L4 (1) 

where, 

L2 = ^Tr{d^Ud^U+) (2) 

is the Skyrme term {U is the chiral field and /vr is interpreted as the pion decay constant), 
and 

L, = ^Trip-'d^U, U+d,U][U+d^'U, U+d^U]) (3) 

is the quadratic term introduced by Skyrme to keep the Skyrmion stable against the 
Derrick instability jZj (e is the Skyrme parameter). Skyrme found that his model contained 
'topologically nontrivial' configurations (extended objects) of the meson fields, namely 
topological solitons, which he identified as baryons. For twenty years the Skyrme model 
was overshadowed by the tremendous success of QCD and only in the early 1980's after 
the establishment of its link to low energy QCD that the model was revived. 

2.1 Skyrmions and QCD 

The success of QCD is limited to the high energy regime, while at low energy it remains 
virtually intractable. The reason for this is that QCD has a running coupling constant ag] 
it is a function of momentum transfer, or distance. At short distances of the order of 0.1 
fm or less (high energy and momentum transfer of several GeV) QCD is characterized by a 
small enough as that it is treated perturbatively. All of the results obtained in this regime 
are consistent with experimental data. This is the phase in which the relevant degrees of 
freedom are quarks and gluons and it is called the asymptotic freedom phase. At large 
distances of the order of 1 fm or more (low energy and momentum transfer of 1 GeV or less) 
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as is of the order of unity and QCD is a nonperturbative theory. This is the confinement 
phase in which quarks are confined inside hadrons and the hadronic degrees of freedom are 
more relevant. This phase, which is the most practical, is the most mathematically com- 
plex. It should provide all properties of hadrons such as masses, sizes, magnetic moments, 
lifetimes, scattering properties and, in principle, all nuclear phenomena. 

The first major step to overcome this problem was taken by 'tHooft [S]. He found that 
in the limit of a large number of colors (large N^), 1/Nc could be used as an expansion 
parameter. In this limit, QCD simplifies a great deal and 'tHooft went on to show that at 
large Nc, QCD is equivalent to a local field theory of mesons and 'glueballs' (bound states 
of gluons, without quarks), with an effective interaction between them of order 1/Nc. The 
second step was taken by Witten jH]. Assuming confinement, he showed that baryons in 
large QCD behave much like solitons in a weakly coupled local field theory of mesons. In 
this limit, baryon masses scale as Nf. = l/g"^, where g is the strength of the meson coupling, 
while baryon sizes are of order 1. Solitons in weakly coupled theories have masses that 
scale as l/g"^ and sizes that tend to constants as g tends to zero. Even though the mesons 
are weakly interacting, the solitons interact strongly as do baryons in QCD. 

The next natural step, it seems, is to derive the effective meson Lagrangian from the 
fundamental QCD Lagrangian. This task, as it turned out, is immensely difficult. Its 
achievement is equivalent to the solution of the intractable original problem of low energy 
nonperturbative QCD in the confinement phase. Nevertheless, the form of the resulting 
effective Lagrangian is being narrowed down under reasonable assumptions. Assuming 
chiral symmetry is spontaneously broken in QCD (with the physical pions as the resulting 
Goldstone bosons, and taking the low energy limit in which one expects the Golstone 
bosons to dominate the dynamics) it has been shown that the first term in the resulting 
low energy effective Lagrangian is |lUj : 

L,^. = ^Tr{d,Ud^U+) + ... (4) 

which is astonishingly similar to the Skyrme term (Eq. |2))! This picture, which emerged 
from large- A^c QCD, is precisely what Skyrme had in mind long before QCD. Further work 
showed that the similarities between the Skyrme model and to what is described as the 
nature of mesons and baryons in large N^. QCD is simply intriguing. 

The Skyrme model, however, as it is build was known to predict an isospin independent 
spin-orbit force with the wrong sign. That is, it predicts a repulsive interaction. 

2.2 Skyrme model and the missing attraction 

The product ansatz for the two-baryon system as suggested in Ref. [Hj (Eq. |21), beyond its 
relative simplicity as compared to other two-baryon field configurations which can be found 
in the literature ^], becomes exact for large N — N separation. Unfortunately, it is 
not the case for the isoscalar component of the spin-orbit force since the standard Skyrme 
model predicts an isospin independent spin-orbit force with the wrong sign. Namely, it 
predicts a repulsive interaction while the phenomenological Bonn potential as the 
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Paris potential gives an attractive one. This came to be known as the problem of the 
missing attraction. Extensions of the Skyrme model consisted on including higher-order 
terms in powers of the derivatives of the pion field ^3 1^ • Expressed in terms of an SU (2) 
matrix U which, as we have said characterizes the pion field, a six-order term corresponding 
to a;-meson exchange [Hj , 

= -^B,{U)B^{U) (5) 

where B^' = e^'^^^^Tr {{dM)U^{daU)U+{d,3U)U+) /2Att^ is the baryon current, the u- 
meson mass and f3^ a dimensionless parameter related to the u tc-/ width, might be a 
good candidate to solve the problem of the N — N isoscalar spin-orbit force. While it was 
believed that the inclusion of such a term leads to the correct sign {attractive interaction) 
for the isoscalar spin-orbit potential [IHl HH], recent calculations (see Ref. [20] for e.g.) 
proved that by taking into account the second part of of the sixth-order term^ the anomaly 
of the Skyrme model {repulsive force instead of attractive) remains. The treatment of the 
spin-orbit part of the two-pion exchange potential within the Skyrme model needed to be 
improved in order to correct the anomaly of that sign. Below we explain how and why the 
dilaton field was suggested as a plausible cure. 



3 The Dilaton field in the Skyrme model 

3.1 The Dilaton contribution and Skyrmion structure 

In Ref. [T^ and [22] the authors explored the idea of coupling the Skyrmion to the dila- 
ton field. This idea to account for a scalar field confines the Skyrmion and provides the 
attractive term missing in the original Skyrme formalism (see OB for more on this). 



3.2 The Skyrmion fluid and Skyrmion stars 

We start by writing the energy of Skyrmions per unit volume (parameterized by the 
density, pv) at finite temperature. In the mean field approximation it is given by (where 
we adopt natural units with h = c = 1), 

^^-^Ep(rip + np) + V;(ao) 

-^£^'"""^^0^0 + ^v^oPv • (6) 

Here, gv is the strength of the coupling of the c<Jo-meson (of mass m^^^; ujq =<u» is 
the mean-field value) to Skyrmions while g]^ represents the isospin degrees of freedom 

•^In Ref. jlH] and the authors considered only one part of the interaction due to the sixth-order 
term in their calculations and omitted the second part which arises from the exchange current |21j . 
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(g'jv = 1 for neutron matter and Qn = 2 for symmetric nuclear matter)^. Ep is the 
contribution to energy of a single Skyrmion while p is the Skyrmion's momentum. Up and 
ftp parameterize the single particle and anti-particle distribution functions. The pressure 
of such an ensemble at T = is then simply given by Py = pyd{Ev / pv) / dpv ■ One finds 
that the contribution of the vector meson field (cuq) to the pressure grows with density 
(tuo cx pv)) and is positive; the dilaton potential Vo-(cro) (where Uo =<o> is the mean-field 
value) gives a negative contribution to the pressure, acting to bind the system into a fluid; 
the SF. Once the EOS describing such a fluid is derived (OB) the corresponding compact 
objects (SSs) are then computed 

SSs in our picture are not boson/soliton stars where the soliton is a global structure over 
the scale of the star but rather form their constituent baryons as topological solitons using 
pions fields. This is fundamentally different from other "Exotic" stars which also follow 
from solutions to an effective non-linear field theory of strong forces (see also [23 121] ; 1211 
and references therein) . In what follows, we conclude this letter by describing astrophysical 
cases plausibly supporting the existence of SSs in nature. 

4 Neutron stars vs Skyrmion stars 

4.1 4U 1636-53 

In Fig.^we compare the M — R relation for Skyrmion stars (OB) to the theoretical M — R 
curve obtained using six recent realistic models for the EOS (UU, BBBl, BBB2, BPAL12, 
Hyp, and K~^). The solid curves labeled SSI and SS2 are for strange stars (the data was 
kindly provided to us by the authors of Ref. [211 )• The triangle depicts the mass-radius 
constraint from fits to X-ray bursts in 4U 1636-53. Inside the triangle is the allowed range 
of M and R which satisfies the compactness constraints as modeled in Ref. j22| (see their 
Figure 4), and clearly favoring stiffer EOSs. Our modeled stars (OB) cross the triangle 
suggestive of 4U 1636-53 as a plausible SS candidate. 

4.2 QPOs and Skyrmion stars: 4U 1820-30 

QPOs are neutron stars emitting X-rays at frequencies of the orbiting accreting matter. 
Such quasi-periodic oscillations (QPO) have been found in 12 binaries of neutron stars 
with low mass companions. If the QPO originate from the innermost stable orbit [211123) 
of the accreting matter, their observed values imply that the accreting neutron star has 
a mass of ~ 2AMq in the case of 4U 1820-30; this would rule out most modern EOSs 
allowing only the stiffest ones. 

^In the energy equation as given by Eq. ((^Jl, the p meson couphng has been omitted which reduces the 
analysis to symmetric nuclear matter only. As such in OB the isospin degrees of freedom ((Jn) was taken 
as a free parameter as to allow for the two regimes ~ pure neutron matter and symmetric matter - to be 
explored. 
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Figure 1: The M — R relation for non- rotating Skyrmion stars (OB) as compared to 
theoretical models of non-rotating neutron stars (UU, BBBl, BBB2, BPAL12, Hyp, and 
K^^) and strange stars (SSI and SS2); see Ref. jlj. The Schwarzschild radius (2GM/c^) is 
shown as a dotted line. Inside the triangle is the allowed range of M and R for 4U 1636-53 
as modeled in Ref. [21! using fits to X-ray bursts. 



SSs is one possibility given that the gravitational mass of the maximum stable non- 
rotating SS has a value of ~ 2.8M0 (OB). For completeness, one should note that even by 
making the modern/recent EOS stiffer at high densities in a smooth way, the maximum 
mass can never exceed 2.3M0 due to the causality condition pUj . 

5 Conclusion 

We gave a brief historical overview of the Skyrme model, its predictions of hadronic inter- 
actions and its interesting connection to QCD. Here we showed how the repulsive term in 
the Skyrme model can be removed by coupling the Skyrmion to the dilaton field. This lead 
to the concept of the Skyrmion fluid and the related hypothetical stars we called Skyrmion 
stars. The stiffness of the Skyrmion fluid allows for SSs to be as massive as 2.8M0. The SSs 
show unique features; for a given mass their radii are in general larger than those of neutron 
stars constructed using modern EOSs. We discussed examples in astrophysics where SSs 
might constitute plausible candidates. Future observations constraining the mass-radius 
plane of compact stars would most likely prove or rule out the existence of SSs in nature. 
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